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ABSTRACT 

We show that the emission from the first, metal- free stars inside Population III ob- 
jects (PopIIIs) are needed to explain the level of fluctuations in the Near Infrared 
Background (NIRB) recently discovered by Kashlinsky et al. (2002), at least at the 
shortest wavelengths. Clustering of (unresolved) Pop Ills can in fact account for the 
entire signal at almost all the ~l-30 arcsec scales probed by observations in the J band. 
Their contribution fades away at shorter frequencies and becomes negligible in the K 
band. "Normal" , highly clustered, (z) ~ 3 galaxies undergoing intense star-formation 
such a those found in the Hubble Deep Fields can 'fill in' this gap and provide for the 
missing signal. It is in fact found that their contribution to the intensity fluctuations 
is the dominant one at A = 2.17/im, while it gradually looses importance in the H 
and J bands. The joint contribution from these two populations of cosmic objects is 
able, within the errors, to reproduce the observed power spectrum in the whole Near 
Infrared range on small angular scales {9 < 200" for Pop III protogalaxies). Signals 
on larger scales detected by other experiments instead require the presence of more 
local sources. 

Key words: galaxies: clustering - galaxies: infrared - cosmology: theory - large-scale 
structure - cosmology: observations 



1 INTRODUCTION 

Different authors have reported the discovery of an excess 
in the Near Infrared Background (NIRB) which cannot be 
accounted for by normal galaxies (see Hauser & Dwek 2001 
and references therein). This unaccounted excess can be well 
fitted by the redshifted light of the first population of stars 
born inside the so-called Population III objects (Pop Ills) 
(Santos, Bromm & Kamionkowski 2001, Salvaterra & Fer- 
rara 2003 [SF03]). More recently, Kashlinsky et al. (2002) 
have observed small-scale NIRB fiuctuations in the J, H 
and K bands, which cannot be due to local or low redshift 
sources. 

The aim of this paper is to show that PopIIIs are indeed 
required to account for the observed small-scale {6 < 200") 
background fluctuations at least in the J band. Their con- 
tribution - dominant for A — 1.25/im - is found to rapidly 
decrease at longer wavelengths until it becomes negligible 
for A = 2.17/im. We will also show that "normal', high 
((2) ~ 3) redshift and highly clustered star-forming galaxies 
of the kind found in the Hubble Deep Fields can provide 
for the 'missing' power in the H and K bands, so that the 
joint contribution from these two populations can reproduce, 
within the errors, the observed power spectrum of intensity 



fluctuations at all NIR wavelengths. Note that, due to the 
high redshifts of the two candidate populations {z ^ 8.8, as 
found by SF03, for PopIIIs and z ^ 2 for "normal" star- 
forming galaxies), this is only true at small angular scales. 
Signals on intermediate-to-large scales detected by other ex- 
periments {9 ~ 0.7°; see e.g. Kashlinsky & Odenwald 2000; 
Matsumoto et al. 2000) are instead expected to require pop- 
ulations of more local sources. 

The outline of the paper is as follows: in Sec. 2 we briefly 
describe the adopted model for the birth and evolution of 
Pop III protogalaxies. Sec. 3 presents our predictions for the 
intensity fluctuations of the NIRB due to the clustering of 
Pop III sources, while Sec. 4 compares such predictions with 
the available data and discusses the results. Finally, Sec. 5 
summarizes our conclusions. 



2 THE MODEL 

SF03 have shown that the NIRB data (see e.g. Hauser & 
Dwek (2001) and references therein) are well fitted by a 
model in which the contribution from PopIIIs is added to 
the "normal" galaxy background light as obtained from the 
Hubble Deep Field (Madau & Pozzetti 2000) or from the 
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Subaru Deep Field (Totani et al. 2001). We now briefly sum- 
marize some of the many aspects of tlie SF03 model. 

The mean specific intensity of the background due to 
PopIIIs I{i>o,zo), as seen at a frequency vo by an observer 
at redshift zo, can be written as 



dz 



(1) 



(Peebles 1993). Here v = + z)/{l + zq), Teff{i^o,zo,z) 
is the effective optical depth at i>o of the IGM between red- 
shift Zo and z, dl/dz is the proper line element, and Zmax 
is the redshift when the sources begin to shine; (z) is the 
comoving specific emissivity in units of erg s~^ Hz~^ cm~^, 
given by 

e,.{z) = l^{z)f^-^ I 
^'■M Jm. 



n{M, z)MdM, 



(2) 



where l^{z) is the specific luminosity of the population at 
redshift z, computed using the spectra for mctal-frcc stars 
obtained by Schaerer (2002); n[M,z) is the comoving num- 
ber density of dark matter lialos of mass M at redshift z 
given by the Press & Schechter formalism (Press & Schechter 
1974). The integral gives the dark matter mass per unit 
volume contained in dark matter hales with mass greater 
than Mmin{z), where Mmin (computed by Fuller & Couch- 
man 2000) is a cutoff mass below which halos cannot form 
stars due to the lack of cooling. and ilb are the total 
matter and baryon density* in units of the critical density 
Pc = 3Hq /StvG. The ratio Qi, /Q,m converts the integral into 
baryonic mass which is then turned into stars with an effi- 
ciency /*. This latter quantity is constrained by the contri- 
bution of the redshifted light of PopIIIs to the NIRB (SF03). 
The PopIII star formation efficiency required to match the 
NIRB data depends essentially only on their initial mass 
function (IMF), being in the range /, = 10% -50% (the top- 
heaviest requiring lowest efficiency). The redshift at which 
the formation of metal-free stars ends is tightly constrained 
to be Zend — 8.8 by the J band data (SF03). 



3 CONTRIBUTION OF CLUSTERING TO 
BACKGROUND FLUCTUATIONS 

The angular correlation function of intensity fluctuations 51 
due to inhomogeneities in the space distribution of unre- 
solved sources (i.e. with fluxes fainter than some threshold 
Sd) is deflned as: 



c{e) = {5i{e' A') 5i{e"A")), 



(3) 



where {6' , 4>') and (6" , (f>") identify two positions on the sky 
separated by an angle 6. The above expression can be writ- 
ten as the sum of two terms, Cp and Cc, the first one due 
to Poisson noise (i.e. fluctuations given by randomly dis- 
tributed objects), and the second one owing to source clus- 
tering (De Zotti et al. 1996 and Magliocchetti ot al. 2001 
for a detailed discussion). In the case of strongly clustered 



* We adopt the 'concordance' model values for the cosmological 
parameters: h = 0.7, Om = 0.3, Oa = 0.7, fl^ = 0.038, ag = 0.9, 
and r = 0.21, where h is the dimensionless Hubble constant. 
Ho = 100ft km s-i Mpc-i 



sources, as the NIR data seems to indicate (Kashlinsky et al. 
2002), the Poisson term Cp is found to be negligible when 
compared to Cc (see e.g. Scott & White 1999). We can 
therefore safely assume C ~ Cc (hereafter simply called 
C{9)), whose expression in the case of Pop III sources is 
given by 



'dx\ 
. dzJ 



du ^g{r,z), 



(4) 



where x = 1{1 + z) is the comoving radial coordinate, r = 
{u^ + x^O^y^^ (for a flat universe and in the small angle 
approximation), ei,(z) the comoving specific emissivity at 
the frequency i^, and zo and Zmax arc defined as in Sec. 2. 

The spatial two-point correlation function of a class of 
gala^xies 'g', ^g{r,z), in general results from a complicated 
interplay between the clustering properties of the underlying 
dark matter and physical processes associated to the forma- 
tion of such galaxies (see e.g. Peacock & Smith 2000; Scocci- 
marro et al. 2001; Magliocchetti & Porciani 2003). However, 
at high enough {z > 0.5 — 1) redshifts, the probability for a 
dark matter halo to host more than one galaxy is negligible 
(see e.g. Somerville et al. 2001), and one can write: 

^s{r, z) = C(n 2)6?ff(M„,i„, z), (5) 

where ^(r, z) is the mass-mass correlation function and 
biff {Mmin, z) the bias associated to all dark matter haloes 
massive enough to host a galaxy at redshift z. 

In equation (5), ^(r, z) has been obtained following the 
work by Peacock & Dodds (1996; see also Moscardini et 
al. 1998 and Magliocchetti et al. 2000), which provides an 
analytical way to derive the trend of the dark matter cor- 
relation function both in the linear and non-linear regime. 
Note that ^(r, z) only depends on the underlying cosmology 
and on the normalization of as- The effective bias factor 
b^ff of all haloes with masses greater than some threshold 
Mmin (quantity which generally may depend on the look- 
back time) is instead obtained by integrating the quantity 
6(M, z) - representing the bias of individual haloes of mass 
M - opportunely weighted by their number density n(M, z): 



bea{z) = 



J^^^JM b{M,z) n{M,z) 
dMn{M,z) ' 



(6) 



where, in this case, the mass function is obtained according 
to the Press & Schechter (1974) approach and we take the 
functional form for b{M, z) from Jing (1998) (but see also 

Scannapieco & Barkana 2002). 

If wo then plug into equation (4) the expressions ob- 
tained by SF03 (see Sec. 2) for the specific emissivity ev{z), 



the halo cutoff mass M„ 



Mjnin{z), and use as the mini- 



mum redshift at which these sources contribute to the back- 
ground the value z^nd ~ 8.8 (corresponding to the epoch 
at which metal-free star formation ends - note that we are 
making the reasonable assumption that no such a source can 
be resolved by current observations), then, by also mating 
use of expressions (5) and (6), we can derive predictions for 
the contribution of the clustering of (unresolved) PopIIIs to 
the background fiuctuations at different wavelengths. 

The cases for A = 1.25pm, A = 1.65//m and A = 2.17/im 
(respectively corresponding to J, H and K bands) are pre- 
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Figure 1. Predictions for the correlation signal owing to inten- 
sity fluctuations of PopIII sources. The solid line represents the 
contribution in the J band, while dashed and dotted lines respec- 
tively illustrate the case for the H and K bands. 

sented in Figure 1 by the solid, dashed and dotted lines. 
Since these wavelengths at the redshifts under exam are al- 
ways greater than the rest- frame Lya, the contribution from 
the e~ f term in equation 4 can be neglected in the cal- 
culation of C{e) (see SF03). 

The first feature to be noticed in the plot is the sharp 
drop of all the curves at ^ ~ 200 arcsec. This is due to 
the fact that such an angular scale corresponds to distances 
r > 8.6 Mpc (with the minimum value corresponding to 
Zend = 8.8 in the adopted cosmology), where the spatial 
correlation function has already steepened from its power- 
law behaviour and rapidly approaches the zero value. As a 
first conclusion of this work we can then say that the clus- 
tering of unresolved PopIII sources cannot account for any 
of the observed fluctuations on scales > 200 arcsec, which 
instead require much more local objects (for the Matsumoto 
et al. 2003 results the maximum acceptable redshift turns 
out to be z ~ 1 — 2). 

The second point to stress is the remarkably different 
amplitude of the intensity fluctuations as evaluated at dif- 
ferent frequencies. More specifically, C{9) is found to de- 
crease by about two orders of magnitude when going from 
A = 1.25 fj,m to A = 2.17 fj,m. As already argued by SF03, 
the reason for this decrement has to be found in the ex- 
tremely strong Lya nebular emission line - responsible for a 
considerable fraction of the PopIII emissivity - which, once 
redshifted to the present time, gives its maximum contri- 
bution in the J band and rapidly disappears at the other 
two wavelengths under exam. We discuss the implications 
of these findings in the next Section. 

4 COMPARISON WITH THE DATA 

Kashlinsky et al. (2002) and Odenwald et al. (2003) have 
recently reported the first detection of small angular scale 




1 10 
k"^ [arcsec] 

Figure 2. Power spectrum of NIRB fluctuations at A = 1.25^m. 
The data are the results obtained by Kashlinsky et al. (2002) 
after atmospheric subtraction; the lower panel reproduces mea- 
surements from their Field 3 (considered by the authors as the 
most representative one), while the upper panel shows the signal 
obtained by averaging all the seven fields. The shaded area corre- 
sponds to uncertainties in the power spectrum arising from both 
statistical errors and fleld-to-fleld variations. In both panels, the 
dashed line represents the predicted contribution to NIRB fluctu- 
ations in the J band as due to clustering of PopIIIs, the dotted line 
that of "normal", highly-clustered, {z) ~ 3 star-forming galaxies 
and the solid line is the sum of both components (see text for 
details). 

fluctuations in the Near Infrared Background. Their mea- 
surements were obtained by using long integration data con- 
structed from 2MASS (Two Micron Sky Survey) observa- 
tions and by then coadding images in order to produce a 
8.6' X 1° field, divided into seven square patches 512" on 
the side. In each patch, individual stars and galaxies were 
removed down to a magnitude limit which slightly varied 
from patch to patch, due to variations in the background 
level and associated noise. The clipping algorithm was ap- 
plied for a mash size of 3 by 3 pixels (corresponding to 9 
square arcsecs). However, in order to test for its stability, 
the above procedure was also repeated with increasing mask 
sizes (up to 5 by 5 pixels) and the corresponding results were 
found in agreement with each other (S. Kashlinsky, private 
communication). Such an approach ensures the effective re- 
moval of entire galaxies at all redshifts but the lowest ones 
in which case more extended objects were identified by eye 
and removed by the clipping algorithm regardless. The re- 
maining diffuse light was then Fourier transformed, with the 
resulting power spectrum showing a positive signal at all 1- 
30 arcsec scales and in all the seven patches. 
The authors estimate contributions to the observed power 
spectrum from possible non-cosmological components such 
as atmospheric fluctuations, remaining Galactic stars, cirrus 
emission, zodiacal light (expected to be the largest source 
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Figure 3. As in Figure 2, but for A = 1.65^m. 
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Figure 4. As in Figure 2, but for A = 2.17/jm. 



for the J, H and K bands. The lower panels show the power 
spectrum as measured at the different wavelengths from 
Field 3, considered by Odenwald et al. (2003) as the most 
representative one due to its depth (limiting magnitudes for 
sources removal: J=20, H=19.2, K=19). The upper panels 
instead report the signal obtained by averaging measure- 
ments in all their seven fields. The shaded areas illustrate 
the levels of uncertainty due to both statistical (Poissonian) 
errors (corresponding to the errorbars in the lower panels) 
and field-to-field variations, found to be the dominant effect 
possibly because of the small area covered by each of the 
patches (cosmic variance). Note that the data has already 
been corrected for the atmospheric contribution. 

In order to compare these measurements with predic- 
tions from Sec. 3, we have converted our results for the inten- 
sity fluctuations C{6) due to PopIIIs into the corresponding 
signal in the angular wavenumber k space, P{k), where 



p{k) = 2tv / c{e) Jo{ke) ede, 



(7) 



with Jo zero-th order Bessel function. 

The resulting power spectrum at the different wave- 
lenghts is illustrated in Figures 2, 3 and 4 by the dashed 
lines. As it is clear by comparing data with model predic- 
tions in both panels of Figure 2, the clustering properties of 
PopIIIs can nicely explain the observed fluctuations in the 
J band, even though there seems to be some discrepancy 
at the smallest angular scales. The situation however gets 
increasingly worse when one goes to longer wavelengths and 
it is found that, while in the H band this class of sources 
can only (possibly) account for the observed spectrum at 
the highest angular scales probed by the data, their contri- 
bution in the K band falls short by at least one order of 
magnitude. Also, especially in the H and K bands, the pre- 
dicted slopes for the power spectrum result too steep when 
compared to the measured ones. 

The missing ingredient to this analysis has to be 
searched amongst "normal galaxies" which can at least par- 
tially account for the extragalactic light emitted at the wave- 
lengths under exam (see Madau & Pozzetti 2000; Totani 
et al. 2001). In their work, Kashlinsky et al. (2002) and 
Odenwald et al. (2003) remove from the maps all resolved 
sources up to K~ 19™ which, they claim, corresponds to 
the removal of possibly all resolvable galaxies to redshifts 
2 ~ 1 (with some allowance made for low-surface brightness 
systems). However, we argue that the population of galaxies 
responsible - together with PopIIIs - for the observed NIRB 
fluctuations must reside at redshifts z > 2 ^. The argument 
supporting this conclusion goes as follows. 

The power spectrum of fluctuations as measured by 
2MASS features a remarkable power-law behaviour down 



of uncertainty in the treatment of NIRB data), instrument 
noise and extinction, and conclude that they all have differ- 
ent slopes and negligible amplitudes when compared to the 
observed signal. Odenwald et al. (2003) and Kashlinsky et 
al. (2002) identify the measured power spectrum as due to 
clustering of faint (K > 18.5-19 mag), high redshift, galaxies 
undergoing significant star-formation. 

Their results (kindly provided by S. Kashlinsky in elec- 
tronic form) are presented in Figures 2, 3 and 4 respectively 



t The possibility for a non-negligible part of the unaccounted sig- 
nal to be due to residual disk emission from low-redshift galaxies 
whose nuclei have been removed can be discarded with a high 
confidence level both on the basis of the arguments introduced in 
this section re the clipping procedure and also by noticing that 
(relatively local) galaxies up to the clipping threshold K~ 19™ 
contribute only a very small fraction of the CIB at this wave- 
length (see e.g. Kashlinsky et al. 2002). So, even if some (outer) 
parts of these galaxies were to remain, their contribution to the 
observed signal would be extremely small. 
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to angular scales of 1 arcsec. Under the hypothesis for this 
signal to be generated by large-scale structure, one has 
that - since the correlation function of dark matter haloes 
(see equation 5) must drop to -1 on distances smaller than 
two virial radii duo to spatial halo-halo exclusion (sec e.g. 
Magliocchetti & Porciani 2003) - the power-law trend ex- 
hibited by the data down to 1 arcsec sets a strong limit on 
the minimum redshift at which these "normal galaxies" ap- 
pear. More in detail, one has that for masses of the order 
of 10*A/q (value taken as a reasonable lower limit for the 
kind of objects under exam) and the cosmology adopted in 
this paper, the condition for a power-law behaviour to hold 
at the minimum scales probed by the data is only met for 
redshifts z >2. 

What are then these z > 2 "normal galaxies", faint 
enough not to be resolved by a K< 19 survey, which can 
significantly contribute to the observed fluctuations of the 
NIRB, especially in the H and K bands? A natural candi- 
date can be found in the population of high-z, star-forming 
galaxies detected in the Hubble Deep Fields to 7ab=29. 
Magliocchetti & Maddox (1999) have measured the cluster- 
ing signal for this class of galaxies and shown that they are 
indeed strongly clustered, with an amplitude A ~ 7 x lO^'* 
(where the angular correlation function was parametrized as 
w{9) = AO'"'^), constant within the errors at all redshifts 
2<«<4.8. 

The contribution of this class of galaxies to intensity 
fluctuations in the NIRB can then be written as Cgai(9) — 
I^w{6), where 7b is the residual background intensity after 
source removal. Suitable values for Ib at the wavelengths 
relevant to this work can be obtained by taking the esti- 
mates for the total extragalactic light emitted in the J, H 
and K bands from Madau & Pozzetti (2000), and by then 
subtracting the contributions from resolved sources as given 
in Kashlinsky et al. (2002). Note that the resulting 7b 's, can 
in principle only be considered as upper limits, since the 
Kashlinsky and collaborators figures were given by assum- 
ing source removal up to « ~ 1, while we are considering 
background fluctuations made by z > 2 galaxies. 

The power spectrum of intensity fluctuations in the 
NIRB obtained for these "normal" galaxies in the J, H and 
K bands is shown in Figures 2, 3 and 4 by the dotted lines. 
The solid lines instead illustrate the total signal steaming 
from the sum of the contribution of this class of sources 
with the one deriving from PopIIIs. "Normal", faint, high- 
redshift galaxies are able to reproduce the observed power 
spectrum at all angular scales in the K band, both in its am- 
plitude and its slope. Their contribution however becomes 
less and less important as one moves towards shorter wave- 
lengths and - at A = 1.25/im - not only can just account 
for ^ 10% of the observed signal, but also has a different 
(shallower) dependence on the angular scale than what the 
data indicates. 

It is then found that, since PopIIIs and "normal" galaxies 
seem to be responsible for the observed fluctuations with 
a different relevance at the different wavelengths (PopIIIs 
mainly appearing in the J band, "normal" galaxies giving 
the strongest contribution in the K band), their sum is able 
to reproduce the data within the errors at all wavelengths 
(with some allowance made in the H band). In more detail, 
by considering both populations we can give account for 
the observed (roughly constant) amplitude of the observed 



power spectrum at the different wavelengths, and can also 
naturally explain the steepening of its slope when moving to 
shorter wavelengths as due to the intervening contribution 
of PopIII sources. Note that the match is particularly good, 
especially if one considers that there are no free parameters 
in the model. Our predictions do not depend on the exact 
functional form for the chosen IMF and star-formation effi- 
ciency /*, as long as their combination is able to reproduce 
the observed NIRB (SF03). Obviously, one of the main un- 
certainties of this work is associated to the determination of 
the slope of the predicted power spectrum. This is due to 
the lack of solid determinations of the two-point correlation 
function, from both an observational and a theoretical point 
of view, at the high redshifts and small masses under exam. 
The fact that our predictions result - within the errors - in 
agreement with the observed trend however gives us some 
confidence that the assumptions made in the course of our 
analysis are (at least) plausible. 

As a last remark wc note that, even if in the above 
analysis we have considered 7b in the different bands as an 
upper limit of the true background intensity, the fact that 
the data at A = 2.17/Lim is perfectly reproduced by the pop- 
ulation of "normal" galaxies suggests this upper limit to be 
very close to the real value for otherwise the curves would 
have fallen below the observations. Under this assumption, 
there is then no need for the emergence of a third, still- 
unknown, population of objects, since PopIIIs, together with 
high- redshift/star- forming "normal" galaxies, can fully ex- 
plain the observed level of fluctuations. 



5 CONCLUSIONS 

In this paper we have used the model of SF03 to derive 
the expected contribution of the clustering of (unresolved) 
PopIIIs to the intensity fluctuations of the NIRB. Compar- 
isons with the recent observations by Odenwald et al. (2003) 
show that this class of sources can indeed account for the 
observed level of fluctuations at almost all the ~l-30 arc- 
sec scales probed by the data, at least in the J band where 
the LyQ nebular emission line - responsible for most of the 
PopIII emissivity - once redshifted to the present time gives 
its maximum contribution. Their relevance however rapidly 
fades away for shorter frequencies and becomes negligible in 
the K band. 

"Normal" , highly clustered, 2 < 2 < 5 galaxies undergo- 
ing intense star-formation such a those found in the Hubble 
Deep Fields can 'fill in' this gap and provide for the missing 
signal. If we in fact assume these objects to produce a clus- 
tering signal as measured by Magliocchetti & Maddox (1999) 
and associate them to a background intensity in in the J, H 
and K bands as taken from Madau & Pozzetti (2000) (af- 
ter source removal), we find that their contribution to the 
intensity fluctuations is the dominant one at A = 2.17/im, 
while it gradually looses importance in the H and J bands. 

The sum of the two PopIIIs and "normal" galaxies com- 
ponents is able, within the errors, to reproduce the observed 
power spectrum of intensity fluctuations in the whole Near 
Infrared range on small angular scales {9 < 200" for PopIII 
sources). The match is particularly good, especially if one 
considers that our model is not associated to any free pa- 
rameter. Signals on larger scales as detected by other ex- 
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periments (see Hauser & Dwek 2001 for review) are instead 
expected to require populations of more local sources. 

The present results strongly support the scenario pre- 
sented by Schneider et al. (2002) (see also Scannapieco, 
Schneider & Ferrara 2003) , in which the first stars whore bi- 
ased towards high masses, thus producing copious amounts 
of ionizing photons and heavy elements which caused cos- 
mic reionization and enrichment of the otherwise pristine 
intergalactic medium. It is only when the gas metallicity 
reached a critical value (still subject to uncertainties due to 
the complicated physics of molecular cloud fragmentation) 
in the range Z = 10~® — 10~*Zq that a relatively abrupt 
transition to the 'normal' star formation mode took place. 
This epoch has left a remarkable imprint both in the in- 
tensity and the fluctuations of the NIRB that we have here 
isolated for the first time. 
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